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Abstract. We bound the location of roots of polynomials that have nonnegative coefficients 
with respect to a fixed but arbitrary basis of the vector space of polynomials of degree at 
most d. For this, we interpret the basis polynomials as vector fields in the real plane, and 
at each point in the plane analyze the combinatorics of the Gale dual vector configuration. 
This approach permits us to incorporate arbitrary linear equations and inequalities among 
the coefficients in a unified manner to obtain more precise bounds on the location of roots. 
We apply our technique to bound the location of roots of Ehrhart and chromatic polynomials. 
Finally, we give an explanation for the clustering seen in plots of roots of random polynomials. 



The Ehrhart polynomial of a d-dimensional lattice polytope Q is a real polynomial of 
degree d, which has the following two representations: 



Here we chose the letter z for the independent variable in order to emphasize that we think 
of iq as a polynomial defined over the complex numbers. The coefficients cq, c^-i and q in 
the first representation are positive, while the others generally can vanish or take on either 
sign. In contrast, a famous theorem of Stanley [11] asserts that all coefficients aj of iq in the 
latter representation are nonnegative, a, > for {) <i < d. 

Such nonnegativity information is also available for other combinatorially defined poly- 
nomials, a case in point being the chromatic polynomial of a graph (cf. Proposition 16.11 
below). An early example of how combinatorial information might be gleaned from studying 
roots of such polynomials is the Birkhojf-Lewis Conjecture, which asserts that no chromatic 
polynomial has a root in the real interval [4, oo). Somewhat ironically, even though it was 
formulated as a new inroad towards settling the Four Color Conjecture (which it implies), the 
latter is now a Theorem, while the former is still open. Nevertheless, since at least 1965 [7], 
the complex roots of chromatic polynomials have received close scrutiny. A well-known recent 
result by Sokal [10] states that their complex roots are dense in the entire complex plane, if 
one allows arbitrarily large graphs. He was motivated by applications in physics to the Potts 
model partition function. 

Coming back to Ehrhart polynomials, first bounds obtained in [I] on the location of the 
roots of iq for fixed d were substantially improved by Braun [3] and Braun & Develin [1]. All 
of these papers use the nonnegativity of the Oj's, but Braun's crucial new insight is to think 
of the value iq{z) at each z G C as a linear combination with nonnegative coefficients of the 
d + 1 complex numbers bi = bi{z) = (^'''^~*)- In particular, for zq to be a zero of iq, there 
must be a nonnegative linear combination of the hi{zo) that sums to zero. 



1. Introduction 
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In this paper, we extend and generalize Braun's bounds on the location of roots for the 
binomial coefficient basis. We propose a unified approach using Gale duality to bound the 
location of roots, that 

• works in exactly the same way for all bases of the vector space of polynomials of 
degree at most d (Theorem I3.3P , and 

• allows one to incorporate arbitrary additional linear equations and inequalities be- 
tween the coefficients aj beyond mere nonnegativity (Theorem 15. ip . This is applied 
in Section [6] to the case of Ehrhart and chromatic polynomials (Figures [6] and [7]) . 

We apply our approach in Section [3] to explicitly bound the location of the roots of poly- 
nomials with nonnegative coefficients with respect to four common bases of P^; the detailed 
treatment of the binomial coefficient basis comprises Section [H Throughout, we focus on 
bounding the location of the non-real roots, as the case of real roots is much more straight- 
forward (Observation 13.4]) . 

In Section [71 we use our method to explain why the roots of "random" polynomials with 
nonnegative coefficients (for a suitable meaning of "random") tend to clump together, by 
tracing this behavior back to properties of the basis polynomials (Figures [9l and [T0|) . 

1.1. Sketch of the method. Let B = {bo, . . . , bd} be any basis of Pd, the ((i-l-l)-dimensional 
vector space of real polynomials of degree at most d in one variable. 

• We regard B as a collection of vector fields: for each complex number z € C, the 
basis elements bo^z), . . . ,bd{z) define a configuration 13{z) = [wo{z), . . . ,Wd{z)) of 
real vectors Wj{z) = (Rebj{z),linbj{z))'^ in the plane R^. This point of view con- 
verts the algebraic problem of bounding the location of roots of a polynomial into a 
combinatorial problem concerning the discrete geometry of vector configurations. 

• We analyze the combinatorics of B{z) in terms of the Gale dual configuration B*{z). 
In particular, there exists a polynomial / = Yli=o with nonnegative coefficients 
Oj > and a root at z = zq whenever the vector configuration B{zq) has a nonnegative 
circuit, and this occurs whenever B*{zq) has a nonnegative cocircuit. 

The important point is that we obtain a semi-explicit expression for B* for any 
basis of Pd, not just the binomial coefficient basis. In fact, for the power basis bi = z*, 
the rising and falling factorial bases bi = , z-, and the binomial coefficient basis 
bi = (^^^~*) we can make the Gale dual completely explicit. 

• In concrete situations one often has more information about /. Gale duality naturally 
allows to incorporate any linear equations and inequalities on the coefficients, and in 
some cases this leads to additional restrictions on the location of roots. 

• As an illustration, we show how the inequality a,^ < oq -I- ai that is valid for Ehrhart 
polynomials further constrains the location of the roots of iq. We also study the case 
of chromatic polynomials, for which Brenti [5] has shown the nonnegativity of the 
coefficients with respect to the binomial coefficient basis. 

• Braun & Develin [4] derive an implicit equation for a curve C bounding the possible 
locations of roots of / = X^j^o ^« (^"""d"') ' method gives an explicit equation 
for a real algebraic curve whose outermost oval is precisely C. 

It is instructive to visualize the vector fields wq, . . . ,Wd for the binomial coefficient basis, 
i.e., when bj{z) = (^'''^~''') = Rj{z) +ilj{z); recall that Wj{z) = (^Rj{z), Ij{z)) . 
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Figure 1. The values of {(^"^^ *) : < i < d} at different points in the 
complex plane, for d = 6. All vectors are normalized to the same length. In 
gray, the locus of points where two vectors become collinear. 

From Figure [H it appears that at points far away from the origin the vectors Wi are all 
"acute", i.e., contained in a half-plane (that varies from point to point), while closer to the 
origin they positively span the entire space. If true in general, this would imply that far away 
from the origin, f cannot have any roots. 

The detailed analysis (and proof) of this observation will take up the bulk of the paper, 
Sections [2] to [6l and in this special case may be summarized as follows: 

Theorem 14.111 Let d be a positive integer and the set of complex, non-real numbers that 
are zeros of non-identically vanishing polynomials of the form 

j=0 ^ ^ 

with aj > for j = 0, . . . ,d. Then is the set of non-real points in the region bounded by the 
outermost oval of the real algebraic curve of degree d — 1 in the complex plane with equation 

z — z ' 

where ~ denotes complex conjugation. This bound is tight, in the sense that any point inside Z^ 
is a root of some such f. Moreover, there is an explicit representation of this equation as the 
determinant of a tridiagonal matrix; see ^ and Proposition \2.4\ 
The real roots of any such f all lie in the real interval [—d, d — 1]. 

From contemplating Figure [U a naive strategy for bounding the locations of the roots 
comes to mind: First, try to prove that for "far away" z the Wi{z) positively span a convex 
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pointed 2-dimensional cone r. Then determine the generators Wk{z), wi{z) of its facets, and 
the locus C of all z G C for which these facet vectors "tip over", i.e., become collinear. By 
continuity, for zq inside C the origin is a nonnegative linear combination of the u^j, and thus 
zq is a possible root. 

The alert reader will perhaps have lost track of even the number of holes in this argument! 
As a sample, it is a priori not clear (but true, at least for the binomial coefficient basis) that 
the Wi{z) in fact span a pointed cone for all z of large enough absolute value. It is even less 
clear (but true in this case) that the vectors spanning facets of r far away from the origin 
will still define facets just before r ceases to be convex closer to the origin. Furthermore, the 
locus C might (and does) have multiple components, suggesting that one has to exercise more 
care when talking about points zq "inside" C. 

However, the real problem with this approach lies with the fact that the locus of collinearity 

of Wi{z) and Wj{z) is the vanishing locus of the determinant Aj^ = ^' , and evaluating 

this polynomial explicitly quickly becomes a daunting task; moreover, it is not at all clear 
how the knowledge of Ajj for any particular basis would help for other bases of Pd- 
We now present our method that overcomes all these obstacles. 



2. Gale duality 

2.1. Overview. Consider a polynomial / = Yli=o^i^i '^^ degree d, expanded with respect to 
a basis B = {bo, . . . , b^} of P^, the vector space of all polynomials in one complex variable of 
degree at most d. For the moment, we will focus on the complex, non-real roots of /. To find 
these, rewrite the real and complex parts of the condition f[z) = in the form 



(1) 



Ri 
h 



Rd 

Id 



( oo\ 

ai 

\adj 



0, 



where Rj = Rj{x, y) and Ij = Ij{x, y) stand for the real and imaginary parts of the polynomial 
bj{x + iy). 

As suggested in the Introduction, we now regard each basis element 6j not as a complex 
polynomial, but as a real vector Wi{x,y) = {Ri,Ii)^ G I^? . Then there exists some polyno- 
mial / with a root aX z = x + iy \i and only if there exist real coefficients oq, . . . , fld with 

d 

'^aiWi{x,y) = 0. 



If we impose the additional restriction that the Oj be nonnegative but not all zero, this is only 
possible if the positive span of the Wi includes the origin. Among all such linear combinations 
summing to zero, we now consider only support-minimal ones, i.e., those with the minimum 
number of nonzero coefficients Oj. In oriented matroid terminology, the ordered collection a 
of signs of the coefficients of such a support-minimal linear combination is called a circuit of 
the (full-dimensional) vector configuration W = [wq, . . . ,Wd) C R^. To proceed, we regard 
W as a.2x {d + l)-matrix. A Gale dual vector configuration W = {wq, . . . ,Wd) C R'^"^ of W 
is the ordered set of rows of any matrix, also called W, whose columns form a basis for the 
(row) kernel of the matrix W, so that WW = [12]. Gale duality states that the collection 
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of signs (T is a cocircuit of W. This means that there exists a linear form g on R'^"^ with 
{sign g{wi) : i = 0, . . . ,d) = a. 

Clearly, any circuit of W has either two or three non-zero entries (unless it is the zero 
circuit, which we exclude from the discussion). Because zq is a root of / if and only if there 
exists a nonnegative circuit of W{zo) = {wqIzq), . . . ,Wd{zo)) C R^, by Gale duality this 
happens if and only if there exists a cocircuit of W{zo) = (Wo{zo), . . . ,Wd{zo)) C IR,"'"^, i.e., if 
and only if there exists a linear form on R"'"^ that vanishes on all of the Wi except for either 
two or three of them, and on those evaluates to the same sign. Geometrically, there must 
exist a linear hyperplane in IR,'^"^ that contains all vectors Wi except for two or three, and has 
those on the same side. 

Thus, we have traded the search for the locus of two collinear vectors among the Wi € 
(a problem involving only two pieces of input data) for the task of finding a Gale dual W 
in the much higher-dimensional space R'^"^, and hyperplanes passing through almost all of 
the Wi — a problem involving almost the entire input! 

That this is not crazy, but instead effective, is explained by the fact that passing to the 
higher-dimensional representation is possible in great generality, and moreover greatly sim- 
plifies the structure of the problem; see Proposition 12.11 below. 

2.2. Implementation. Let B = {bi : < i < d} he any basis of P^- 

2.2.1. The Gale dual. Form the matrix 



where Rj = Rj{x,y) and Ij = Ij{x,y) denote the real and imaginary part of the complex 
polynomial bj = bj(x + iy). The rank of W is 2, so any Gale dual matrix to has size 
{d + 1) X (d — 1). The following proposition gives an explicit representative for W involving 
polynomials Pk,Qk,rk that depend on the basis B. For four especially relevant bases, we will 
make the Gale dual W completely explicit. These bases are: 

• The power basis, where bi = z^; 

• the falling factorial basis, where bi = z- = z{z — 1) ■ ■ ■ {z — i + 1); 

• the rising factorial basis, where bi = z^ = z{z -|- 1) • • • (z -|- i — 1); and 

• the binomial coefficient basis, where bi = {^~^d~^)' 



Proposition 2.1. A Gale dual matrix to W may be chosen to have exactly three non-zero 
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Moreover, its entries may be chosen to lie in M.[x,y]. For the four bases considered, we may 
choose the following explicit values: 



z 
z'- 



Pk 



Note that in the last row, 



x"^ + y"^ 
{x — + 
(x + + y2 
(x — + y"^ 

= 2(x-(i 



Qk 



2x 

2{x -k)-l 
2{x + k) + l 
Pk + rk- d{d - 1) 
2 . - ,2 (f^l 



rk 



1 
1 
1 

Pk+l~d 



^)r + 2y' 



Proof. We first prove that the matrix W can be chosen to have the displayed triple band 
structure regardless of the basis B chosen for P^- For this, define the rational functions 

1; specific values for gk become apparent from the relations 



5^ = ^ G ]R,(z) for < A; < d 



yk + l 



z- '-=z- z\ zttl ={z- k)z^, z^+'^ = {z + k)z^ and {'^'^'J"^) = 
The triple {pk^Qki^k) lists nontrivial coefficients of a real syzygy 



Pkbk + Qkbk+i + rkbk+2 = bk {pk + gkQk + Qkgk+iTk) = 



whenever 



1 Regk Res-fc^fc+A 
\mgk Imgkgk+iJ 



But the displayed matrix with entries in lR,(x,y), call it M, obviously has rank at least 1, 
and rank 2 whenever Im gk{x + iy) ^ 0, so that such triples certainly exist. Moreover, by 
multiplying with a common denominator we may assume Pk,Qkjfk € lR,[x,y], and so the 
relations pkbk + Qkbk+i + rkbk+2 = imply that W is in fact a Gale dual of W. The concrete 
syzygies listed above arise by choosing explicit bases for ker M. □ 

Remark 2.2. Another interesting case is that of polynomials with symmetric coefficients. 



For instance, if / = Yli 



. a,; 



*) and 



a,; 



ad- 



we may expand / in the basis B 



{ {^~^d~^) ~'~ (^d *) • — ^ — |_f J } of the vector space of polynomials with symmetric coefficients 
in the binomial coefficient basis. However, the coefficients of syzygies of these bk do not appear 
to be as simple as the ones listed in Proposition [2TTJ For example, a typical coefficient (namely, 
qi for d = 8) reads 



-8(((x + aif + + aa)' + y^) + 71) (((^ + Pi? + + 1^2? + y^) + 72) + 7 

where ai, 02 are the roots oi — a + pi = (so that ai + as 
P"^ - P + P2 = 0, pi,p2 are the roots of - - ^ = 0, 71 + 72 = 
and 7 = We will not pursue this basis further in this paper. 



1), Pi,P2 are the roots of 

_ 135 ^, _ 135/1 61 



649' 



71 



2.2.2. The determinants. Recall that two vectors u!j{z), Wk{z) become collinear at some point 
2; G C whenever there exists a circuit of the vector configuration VF(2:) with exactly two non- 
zero entries. By Gale duality, this means that the Gale dual vector configuration VF(2:) has 
a cocircuit with support 2, i.e., the determinant of the matrix obtained by deleting two rows 
from W vanishes. Our approach rests on the fact that we can give fairly explicit expressions 
for these determinants for the four bases considered here. 
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Lemma 2.3. Let W^j^f^-^ = W (^j^j^-j{x,y) denote the square matrix obtained by deleting rows j 
and k from W , where < j < k < d (so that we number the rows from to d). Then 

where Dj k = Dj j.[x,y) is the determinant of the tridiagonal matrix 



(-Ij 



(3) 



V 



rk-3 



\ 



Pk~2 

-qk-2l 



Here -Djj+i := 1, and the leading resp. trailing products are 1 if j 
particular, Djj^2 = —Qj- 



resp. k = d. In 



Proof. The matrix fc) decomposes into three blocks, whose determinants yield the stated 
formula, and two additional elements rj_i and Pk~i that do not contribute to det VF(j^fc). □ 

Proposition 2.4. Set z = x + iy and z = x — iy. Then 



Dj,k{x,y) 



■{fj,k{z) - fj,k{z)), 



where the polynomials fj^k{z) are given in the following table: 



bi 




z' 


zk-3 




z^ 


{z-k + l)--- 


{z - j) 




{z + k-l)--- 


(z+j) 




\{z-k + l)---{z-j){z + d 


-k + l)---{z + d-j) 



The Dj^k are real polynomials with even degrees in y. 

Proof. It is well known that the determinant D„ of an n x n tridiagonal matrix A = (aij) 
satisfies the three-term recursion relation Dn = OnnDn-i — an,n-ian-i,nDn-2- Solving this re- 
cursion for the matrix from Lemma [2. 31 with the values from Proposition 12.11 and the boundary 
conditions Dj,j+i = 1 and -Dj,j+2 = —<lj yields the stated expressions. □ 

2.3. The real case. Up to now, we have only considered complex, non-real roots of /. The 
case of real roots is much simpler, and the machinery used for complex roots specializes in a 
straightforward way to the real case. If we regard both / and the 6j as polynomials in one 
real variable, the matrix W = reduces to the single row = (6oi ■ ■ ■ ,b(i). 

Proposition 2.5. A basis for the kernel ofW^ is given by the columns of the matrix 



W 
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of size {d+ 1) X d. The determinant of the matrix obtained by deleting row j from W is 

det Wj = {-lybjU for j = 0,l,...,d, 
with n = 6162 • • • bd-i ■ □ 

3. Bounding the location of roots 

We first treat the case of complex, non-real roots. For each ordered triple of indices i,j, k 
with 0<i<j<k<d, denote by Hiji. the hyperplane in IR''"^ spanned by the rows of the 
matrix VF^jj- ^), obtained by deleting the rows Wi,Wj,Wk from W. 

Definition 3.1. Sij^k is the set of all 2; = x + «?/ G C such that Hij ^ = Hij ^ix, y) induces a 
nonnegative cocircuit, i.e., the vectors Wi = Wi{x, y),Wj = Wj{x, y),Wk = Wk{x, y) all (weakly) 
lie on the same side of Hij^k{x,y)- 

The sets Sij^k are crucial for our purposes for the following reason: If 2; € <Sij,kj then 
the corresponding Gale primal vectors Wi,Wj,Wk form a nonnegative circuit, and thus yield 
a nonnegative combination of all w^s that sums to zero; in other words, there exists some 
polynomial / with nonnegative coefficients in the chosen basis B that has a zero at z. On the 
other hand, if 2; ^ <Sij,k, we can only conclude that the three particular Gale primal vectors 
Wi,Wj,Wk do not form a circuit, and so are not responsible for the possible zero z of /. 

Proposition 3.2. For 0<i<j<k<d and z G C \ R, let Cij^kiz) be the set of signs 

|(-l)*signdetW(j-fc)(x,y), (-l)-'"^ sign det W(,,fc)(x, y), signdetW^,j){x,y)^ 

= |(-l)*signl)j-fe(x,y), (-1)^-^ sign A,fc(a;, y), (-1)'="^ sign A,i(2;, y)} . 
Then each Sij^k C is a semialgebraic set defined as the locus of all (x,y) such that 

{±1} (Tij^kix + iy). 

Proof. We obtain a linear form ipi,j^k on R*^"^ whose vanishing locus is the hyperplane Hij^k 
by adding a first row of variables xi, . . . ,Xrf_i to W(^ij^k) expanding the determinant 
of that square matrix along the first row. The value of <Pij^k on Wi, say, is given by the 
sign (—1)* of the permutation that interchanges rows and i in the matrix W^j^i^y times 
detW^^k)- For Hi,j,k to define a (positive or negative) cocircuit, the signs obtained in this 
way for Wi, wj and Wk must agree. Finally, by Lemma [2T3l the signs of det W^^k) and Dj^k 
agree except perhaps on the real axis (on the vanishing locus of the pkS and r^'s), and we 
may assume that j — i >2 and k — j >2. □ 

In summary: 

Theorem 3.3. Let f be a polynomial of degree d with nonnegative coefficients with respect 
to some basis of the vector space Pd- Then the set of non-real roots of f is contained in the 
union of the semialgebraic sets Sij^k, for 0<i<j<k<d. Put differently, if 

{-1,1} C {{-iyD,M.zo), {-iy+'D,M.zo), (-l)'=A,,(^o)} 
for each triple {i,j, k) with 0<i<j<k<d, then zq is not a root of f . □ 

After a short discussion of the real case, we will apply this result to our four representative 
bases. We only discuss the power basis and binomial coefficient basis in any detail, as the 
procedure for the rising and falling factorial bases is almost exactly the same. 
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3.1. The real case. A real number x G IR, is a real root of some polynomial / with non- 
negative coefficients with respect to a fixed basis if and only if either 2; is a root of some basis 
polynomial, or two basis polynomials differ in sign when evaluated at x. In other words: 

Observation 3.4. Let f be a polynomial of degree d with nonnegative coefficients with respect 
to some basis {bo, . . . ,bd} of P^. Then the locus of possible real roots of f is the set of x G IR, 
for which bi(x)bj{x) < for some i ^ j- □ 

For the sake of completeness, and in response to the query of one of the referees, we briefly 
rederive this result using our framework of Gale transforms. 

Proof. In complete analogy to the complex case, denote for 0<i<j<dhy W^(ij) the 
matrix obtained by deleting rows Wi and Wj from W^, by Hij the hyperplane in M!^ spanned 
by the rows of W^(jj), and by Sij C ]R the set of all x G IR, such that Hij{x) induces a positive 
cocircuit, i.e., the vectors Wi = Wi{x) and Wj = Wj{x) lie on the same side of Hij. 

To find these cocircuits explicitly, build a linear form ipij on IR'^ that defines Hij by adding a 

first row {xi, . . . , Xd) of variables to W^j and expanding the determinant of that square matrix 

along the first row. Just as in the proof of Proposition 13.21 (fij{wi) = (— l)*detM^j = ^jll 

and ipij{wj) = (— 1)-'"^ det = —bjH by Proposition 12.51 In consequence, Sij is the locus 
of points X G IR such that and bj{x) differ in sign. This finishes the proof. □ 

3.2. The pov^er basis. For bi = z^, we set n = k — j — 1, write -D„ for -Dj,fc, and substitute 
z = re*^ into 

n+l _ ^n+l 2i(n+l)e _ 1 

Dn = i-ir ^ = (-l)V-e— ^ ^.^ . 

z — z e^*" — 1 

This vanishes iff = 7r//(n + 1) for integer I with 1 < / < 2n + 1 and / 7^ n + l. The zero 
locus of Dn thus consists of n lines through the origin, the ones closest to the x-axis having 
angles 9 = We conclude that Dn has the same sign throughout the entire open sector 

Zn+i = {z G C : — < argz < ^^}. By substituting a positive, real value of z into 
Dn = (-1)" Ei=o -2^'^""^ we determine this sign to be (-1)" = 

For z Zd and 0<i<j<k<d the set of signs of the polynomials in Proposition 13.21 is 

= {{-iy+^+''+\ (-l)*+^■+^ {-iy+^+'^-^} = {±1}. 

This implies Sij^k = 0) and thus Theorem 13.31 recovers the classical result that a polynomial 
of degree d with positive coefficients in the power basis has no zeros in Z^,; of course, this 
includes the case of real roots. 

3.3. Rising and falling factorial basis. In both cases, the polynomials fj^k from Propo- 
sition [23] have the form = 11^=1^ (-^ ~ ^i)^ with n = k — j — 1 and Oj = j — 1 + i for 
the falling powers and Oi = —{j — 1 + i) in the case of the rising powers. The transform 
2; I— > z lb '^'^2^^ remedies this asymmetry, where we choose the ' — ' sign for 6j = z- and the 

'+' sign for bi = z"^ . The then become integers or half- integers in the range it^. 

Using the same type of analysis as will be detailed in Section U] for the binomial coefficient 
basis, one can prove that the zero locus Vj^ = {2; G C : Dj^^i^) = 0} is smooth everywhere, 
that one component intersects the real axis between each pair of adjacent a^'s, and that far 
away from the origin "Dj^k approaches the arrangement of lines through the origin with slopes 
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±^^, . . . , ±^^; cf. Figure [21 We will not enter into the details here, but instead treat the 
remaining basis in a separate section. 




Figure 2. The locus Vq iq in the case of the rising factorial basis. 



4. The binomial coefficient basis 

We first get the real case out of the way: The basis polynomials all have the same sign 
outside the closed interval [—d,d — 1], and at each point inside this interval there are two 
basis polynomials that evaluate to opposite signs. By Observation 13.41 [—d,d— 1] is exactly 
the set of possible real roots. 

For the non-real roots, as before we pass to an adapted coordinate system with respect to 
which the vanishing locus of Dj^^ is centro-symmetric, by replacing 

(4) z ^ z' + {k + j -d-l)/2 

in d ■ fj^kiz)- Writing again z for z' yields 



fe-i 



k+j-d-l-2i 



k-i 



n 



k+j + d-l-2i 



1=] - ' i=j 

Next, we replace i hy i + j in the first product and hy k — 1 — i in the second, to obtain 

k-j 



i=l 



dfj,k{z) = Y{ \ z -i- — ] { z + i + 



A 



where A = d— 1 — k + j. Introducing aj = z + y we obtain 

\n+l "+1 



(5) D,^k{z) = D^{z 



(-ir 



d{z — z) 



W{z- ai){z + ai)-^{z- ai){z + 



where we have set n = k — j — 1 (so that A = d — 2 — n), in accordance with the fact that 
the degree of Dj^k{z) in z is n. 

Before examining the zero locus of Dj^]^[z), we pause to calculate the leading coefficient. 
This result will be used in Section 16.11 
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Lemma 4.1. The leading coefficient of Dj k[z) is 

(6) [2n]Z),-fc(z) = {zz)''-^-\-l)'^-^-Hk-j) = r2"(-l)"(n + l), 

where z = re**^. It is invariant under substitutions of the form z i—>- z + zq, and the sign of 
Dn{z) outside the outermost component ofT>n is (—1)", and +1 inside the innermost one. 

Proof. See the Appendix. □ 

We now treat the zero locus of Dn. First, whenever Dn{z) = 0, 

n+l n+1 n+1 n+1 

arg(z - Oj) + ^ arg(z + Oj) - ^ arg(z - Oi) - ^ arg(z + a^) = 2/7r 

j=l i=l i=l j=l 

for some integer I. Because arg(z it Oj) = — arg(z it Oj), this relation reads 

n+1 n+1 n+1 

(7) ^arg(2; - aj) - ^arg(2; + Oj) = ^ = /vr, 

i=l i=l i=l 

where Oj is the angle under which the segment [—ai, ai] appears as seen from z (cf. Figure [3]). 




Figure 3. The segments [— aj,aj] as seen from z 

We may assume without loss of generality that z lies in the upper half plane, and therefore 
that arg(z — Cj) > arg(2; + Oj) > 0, which implies I > 1. On the other hand, the maximal 
value (n + 1)tt of ([7]) is achieved for real z between — ai and ai, so that / < n for non-real z. 
Prom this, we can draw several conclusions, which we detail in Section 14. 1[ The reader may 
want to just skim this material, and otherwise skip ahead to Section 14.21 where we apply it 
to conclude that the root locus is bounded. 
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4.1. Limiting behavior and global geometric properties of Dj f^. 

Proposition 4.2. When d becomes large with respect to n, the zero locus of Dj^k approaches 
a union of circles passing through zizd/2 and symmetric about the imaginary axis. For I = 
1, . . . ,k — j — 1, these circles have center 

d + 1 — k — j d In 
zi = i - cot 



2 k-j 
and radius 



= ; — • 

2sin-j^ 

Proof. For d large with respect to n, the points itaj fuse to ±a = so that ([7]) reads 

l-K 



a := arg(z — a) — aTg{z + a) 



a = I and equation ([7]) says a = ^. In that limit, Pi thus approaches the circumference with 
center and radius i . For smaller values of d, directly evaluating equation yields 



n + 1 

By elementary geometry, the locus of these points is a union of two circular arcs with the 
specified equations. □ 

Example 4.3. For n = 1, we obtain ai = and 02 = For large d, they approach 

2 and equation ^, — zl in +1,0+ i^, 

and radius 

Di{z) = zz — aia2, 
which describes a circumference of center and radius ^/aia2 = \-\/ d"^ — 1 < ^. 

Proposition 4.4. The plane algebraic curve "Dj k with equation Dj^k{z) = is smooth. The 
only points where it has horizontal tangent vectors lie on the y-axis. 

This is proved in the Appendix. 

Proposition 4.5. All algebraic curves Vj^k consist of n = k—j — 1 nested ovals. The i-th oval 
intersects the real axis inside the union of open intervals ±(0^, aj+i), for i = 1, . . . ,n. 

Proof. Let € 5^ \ 5'^ be a non-real unit vector and p the ray through the origin and cp. At 
each point p of Vj^k H p, the angle sum "^^^i cti takes on some value In among the discrete set 
{tt, . . . , nn}, and therefore this value remains constant on the entire connected component to 
which p belongs. The argument extends to the real axis by smoothness of 

For the second statement, observe that the value of aj = arg(z — aj) — arg(z + aj) increases 
by almost vr as z travels from aj + e + iS to aj — e + i6, for < 6 <^ e <^ 1. □ 

Example 4.6. For 1^0,10 = Dg and d = 10, we obtain the picture of Figure [H 

To continue, we introduce some useful notation. By ([5]), the formula for Z)j ^(z) involves 
the points aj^k;i = ^ + ^{d — 1 — k + j) for l<i<k — j,so that 

f d k — j — 1 d k — j — 1 

(Oj,fc;l) • • • J 0'j,k;k-j) = 2 ' ' ' ' ' 2 2 

We write a{±aj^k;i] z) for the angle under which z G C sees the segment [ — aj^k;i: Oij,k;i 

1, and 

A(j,k;z) = a(±aj,fc;i; -2) for the corresponding angle sum. Moreover, let 

T^j,k;l = {{x,y) eU"^ : A{j,k;x + iy) = In} ior I = I, . . . , k - j - I 
be the l-th oval of Pj,fc) and clT>j^k-i the closure of the region in 

1R2 bounded by Vj^k-i. 
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Figure 4. Left: The curve 'Dq^iq for d = 10; Right: the curves "Dij with 
< i < j < 5 

Remark 4.7. The arrangement of ovals {'Dj^k-i :0<j<k<d,l<l<k — j} has several 
interesting combinatorial properties, which we will not pursue in this paper. Here we would 
only like to point out the triple points where components of X^j,r, I'r.fc and T^j,k intersect. 

Proposition 4.8. Let 0<j<j'<k'<k<d and l<l<k— j — I be integers. 

Then 1^j',k' C clT>j^k-i \ clT>j^k;k-j~i- In particular, all components of all curves 'Dj^k o,Te- 
contained in the topological closure of cWq ^-,! \ clPo,d;d-i- Moreover, for all integers 81,82 
with < (5i < j, < 82 <: d — k and (where appropriate) 1 + 81 + 82 < I < d — 81 — 82, 

(8) cl'Dj_Si,k+S2;l+Si+S2 ^ ^^I^j,k;l - ^1 2?j+5i,A:-52;Z-5i-52 ' 

(9) c\Vj+Si,k^52\l ^ ^^1^j,kd ^ clVj^Si,k+52;l- 

Proof. We first show that T>ji^k'\i C clP^^fc;! for all / with 1 < I < k' — j' — 1. The first 

set consists of all points z € <C such that X]i=i'' 0({^(^j',k';i'i z) = In in the centro-symmetric 
coordinates. Undoing the coordinate change ([4]) yields 

'Df,k';l = {ze(D:a{±{d-k' + l);z)+--- + a{±{d-j'y,z)=lTr}, 

dVj^k;i = {z e<C : a{±{d- k + l);z) -\ \- a{±{d - j); z) > tt}. 

Now the required inclusion is clear, because the first set of points of which the viewing angle 
is taken is a subset of the second one. It remains to prove that T>j^k;k-j-i C cWji^k'-^i for all 
\ <l < k' — j' — 1; proving the extremal case I = k' — j' — l\s sufficient. Thus, we are required 

to show that ^) = (A:-j-l)7r implies Yfmid-k'+i^i'^'^'^^) - (^'-j'-l)^- 

But this is true because the first sum has k — j summands, the second k' — j' summands, and 
removing each of the (A; — j) — (k' — j') = (k — k') + {f — j) > pairs of points from the points 
corresponding to the first summand decreases the total viewing angle by at most tt. 

Similarly, the first inclusion of ([8]) follows because Ylm=d-k+i~52 ^{^''^i z) > {I + 8i + 82)1^ 
implies Ylt^dd-k+i '^(^''^'i ^) — removing 81 + 82 pairs of points, and the second one 

from an appropriate change of variables. Relations ([9]) are proved in exactly the same way. □ 
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Corollary 4.9. For z ^ clPo,d;i; facets of the cone t{z) = r>o(wo(2), • • • , Wd{z)) are the 
rays spanned by wo{z) and Wd{z), and 'Wi{z), . . . ,Wd-i{z) appear in cyclic order inside t(z). 

Proof. The argument of Wj{z) = {^~^d~'') ( SiLi arg {z + i — j)) mod 27r, so that the differ- 
ence of the arguments of Wj{z) and Wj+i{z) equals (3j := arg(2; + d — j) — arg(z — j) mod 27r; 
cf. [4j. If we choose z to have the form z = N + ie, with ^ e > 0, it is not necessary 
to reduce Pj modulo 2tt, and < /5o < /?i < • • • < /S^; we may even achieve Pd < f ) so 
that the total angle subtended by the Wi{z) is strictly less than vr, and wq{z) and Wd{z) span 
the facets of t{z). Now note that two vectors Wi{z),Wj{z) become collinear iff there is a (not 
necessarily positive or negative) circuit involving the two, iff there is such a cocircuit involving 
Wi{z),Wj{z), iff T>j^k{z) = 0. An invocation of Proposition 14.81 finishes the proof. □ 

We close with a lemma regarding the relative orientations of Wj^Wk on T^j k. 

Lemma 4.10. Let z € 'Dj,k\l, md regard Wi{z) = (^'^^~*) ols a vector in R^. Then Wj{z) and 
Wk{z) point in the same direction iff I is even, and in opposite directions iff I is odd: 

sign {wj{z) ■ Wk{z)) = (-1)' for z G D^- fc.;. 

4.2. Conclusion: the root locus is bounded. 

Theorem 4.11. Let f = X]j=o (^'''d""') ^ polynomial of degree d with nonnegative co- 
efficients Oj > with respect to the binomial coefficient basis. Then all non-real roots of f 
are contained in the region clPo,d;i bounded by the outermost oval of the algebraic curve with 
equation DQ^d{z) = 0, and any point inside clPo.d;! arises as a root of some such f . 
The real roots of f all lie in the real interval [—d,d — 1]. 

Proof. The first statement can be proved by a short calculation involving Lemma 14.11 and 
the general tool of Theorem 13. 3i However, we have accumulated enough information about 
the special curves "Dj^^ arising for the binomial coefficient basis to give a direct proof: By 
Corollarv l4.91 the vectors wq = wq{z), . . . ,ujd = Wd{z) are positively spanning for z ^ clPo,(i;i- 

Next, suppose that z € clPo,d;i falls inside the region Sk ■= clPo.fc;! \ cH>o,fc-i;i for some 
k € ¥\ with 2 < k < d. Such a k exists, because cH>o,fc-i;i C clPo,fe;i by ([9]), and cIPq,!;! = 0. 
We claim that in this situation, the vectors wq, Wk-i and Wk are positively spanning. Indeed, 
the locus of points in the complex plane where the combinatorics of this subconfiguration 
changes is exactly X'o.fe-i U I'o.fc) because Cfc.fc-i = 0- Moreover, Po.fc;/ C clPo,fc-i;i for Z > 2 
by dSI), so the boundary of the region Sk is Vq ^-i-i U Vq ^-,!, and the property of the three 
vectors being spanning or not remains constant inside Sk- Since outside of clPo,fc;ij these 
vectors are not positively spanning by Corollarv 14.91 but this changes when crossing dSk, the 
second statement follows. 

Finally, the case of real roots was dealt with at the beginning of the present Section HI □ 

Example 4.12. Let d = 3. Then 





Po 


\ 




-qo 


Pi 




ro 


-qi 


v 





n ) 



qo = 2(x+l)2+2y2_4, = 2x'^ + 2y'^ - 4, po = x^ + y^ pi = (x-lf+y'^, ro = {x + 2f+y^ 
and n = (x + 1)2 +y2. Furthermore, L>o,2 = -qon, -Do,3 = qoqi - Piro, -Di,3 = -poqi, and 
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Dij > otherwise. Now 

Sou = {z : i?i,2 > 0, Z)o,2, Do^i > O}, 5oi3 = {z : ^1,3, 1^0,3, -Do,i < O}, 

5o23 = {z : D2,3 > 0, -Do,3, -1^0,2}, 5i23 = {z : -D2,3 < 0, -2^1,3, -Di,2 < 0}, 

so that by Figure [5] and Theorem 14.111 ah non-real roots of polynomials of degree 3 with 
nonnegative coefficients in the binomial coefficient basis lie in the union of these regions. 




Figure 5. From left to right, the semialgebraic sets 5oi2, 5013, 5023, 5123 
(shaded). Their union equals the entire interior of the bounding curve C, 
which by Theorem 14.111 is precisely the locus of possible non-real roots. 



5. Incorporating additional linear constraints 



5.1. Linear inequalities. Suppose we not only know that the coefficients of a polynomial 
/ = Yli=o ^i^i with respect to some basis B = {bi : i = 0, . . . , d} are nonnegative, but also 
that they satisfy a linear inequality Yli=o'^i'^i — '— case is of course accounted for 

by reversing the signs of the Aj. We use a slack variable s > to rewrite our inequality as 

d 



1=0 



Xitti + s 



0. 





To incorporate this into our Gale dual matrices W and W, we introduce the vector o 
(oo, . . . , a^, s)^ . The analogue Wa = of ([1]) is 

/ao\ 

ai 

and we name the columns of this new W by wq, • • • , Wd+i- We obtain a Gale dual W = W{z) 
of W by appending the row vector 

Wd+i = {ujQ,...,iOd-2) = {-XiPi + Xi+iQi- Xi+2ri ■■ 0<i<d-2) 

to the matrix W from ([2]). For the polynomial f{z) = X]f=o ^«^«(-^) with a,, > and 
Yli=o ^i^i ^ to have a zero at z = zq, the vector d must lie in the column space of W{zq) 
and have nonnegative entries; equivalently, there must exist a vector /x = {fiQ, ■ ■ ■ , A'd-2)"^ with 
W(zo)iJ, = d. Geometrically, we think of /i as the normal vector of a linear hyperplane that 
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leaves all vectors uJj (weakly) on one side. In particular, if the linear inequality is strict (so 
that s > 0), then we are only interested in linear hyperplanes that do not contain Wd+i- 

In general, m > 1 independent linear inequalities yield a {d + 1 + m) x (d — l)-matrix 
W. Consider the configuration of d + 1 + m vectors in IR*^"^ spanned by the rows of W. 
Each {d — 2)-tuple of vectors among these spans a linear hyperplane, and we would like to 
know when the m + 3 remaining vectors all lie on the same side of it. As before, we treat 
strict inequalities by only considering those linear hyperplanes that do not contain any of the 
m "new" vectors Wj, and to simplify the discussion we will focus on these. 

We thus fix an ordered subset J = {ji , . . . , jm+s } = {ji < J2 < is} U {d + 1, . . . , d + m} 
of {0, . . . ,d + m}; this set will index the rows of W not on a linear hyperplane. Next, we 
calculate a linear form ipj whose vanishing locus is the hyperplane spanned by the d — 2 
vectors not indexed by J: it is the determinant of the matrix obtained by deleting from W 
all rows indexed by J, and adding a first row of variables. The sign aj^{z) of ^j{wj-) at a 
point z € C is then obtained by plugging the coordinates of Wj^ = Wj.{z) into these variables, 
i.e., by not deleting the row with index ji, but instead permuting it to the first row and then 
taking the sign of the determinant of the resulting matrix. More precisely, if we denote by Wk 
the matrix obtained from W by deleting the rows indexed by i^' C {0, . . . ,d + m}, then 

(10) aj/z) = i-iy'+'+^ signdet Wj^{j^y{z), for f = 1, . . . , m + 3. 

Writing a{J, z) = {aji{z), . . . , crj ^^^{z)}, we can summarize our discussion as follows: 

Theorem 5.1. Assume that the coefficients of f satisfy m > 1 strict linear inequalities, 
indexed from d + 1 to d + m. Let 

S{J) = {z G C : cj(J,z) = {-1,0} or a(J,z) = {0,+l}}. 

Then the set of roots of f is contained in the union \JjS{J), where J runs through all sets 
of the form {ji, j2, Js} U {d + 1, . . . , d + m} with < ji < 32 < js < d; put differently, if 
{— 1, 1} C (j( J, z) for each such J, then zq is not a root of f . □ 

In the case m = 1 and J = {j, k,l,d + 1}, we obtain from (jlOp that 

(11) aiJ,z) = S^{-iy sign Dk,i, (-1)^+^ sign Z),- (-1)' sign D,- (-l)'^signdet 
Expanding the last determinant along its last row yields 

d-2 

(12) detW{,^k,i} = {-l)''Y.{-iru^i[W]{j,k,i};c, 

c=0 

where [VF]{j,fc,/};c stands for the minor of W obtained by deleting rows j,k,l and column c. 
This formula can be evaluated as follows: 

Lemma 5.2. Let m = 1, 0<j<k<l<d, and < c < d — 2. Then 

_ ipo- ■ -pj-iDj^c+iPc+i- ■ ■Pk-iDk,iri-i- ■ ■rd-2 if < c < k - 1, 
{3,k,l},c . . . pj^^Dj^kn-i ■ ■ ■ '■c-i^c+i./n-i • • • rd-2 if k-l<c<d-2. 

Here we follow the convention that Pa ■ ■ ■ Pb = fa ' ■ ■ = ^ if a > b, but Da^b = for a > b. 
In particular, = for < c < j — 1 and I — l<c<d — 2. 
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Proof. In each case, W k-c decomposes into square blocks on the diagonal whose determinants 
yield the stated expressions. The elements outside these blocks do not contribute to [W^]x;c5 
because the determinant of a block matrix of the form ( ^ {5 ) or ( ^ ^ ) is det A det D. □ 

To recapitulate, additional linear inequalities can only restrict further the location of pos- 
sible roots of /. If the vectors Wi{zQ).,Wj{zQ),Wk{zQ) € S do not witness a possible root 
of/, in other words {±1} C [{-ly s\gxi.Dk,i{zo), (-1)''+^ signD^- ((zq), (-1)' signi:)j,fc(zo)}, 
nothing changes after incorporating the additional sign (—1)'^ sign det W^j k i^{zQ): the vectors 
Wi{zo).,Wj{zQ),Wk{zQ),wi{zQ) do still not witness a root of / at zq. If, on the other hand, the 
new sign is different from the old ones, there is "one reason less" for zq to be a root. 

5.2. Linear equations. If the coefficients of / satisfy m independent linear equations of 
the form Yl'l=Q\ai = (corresponding to the case s = 0), the d — \ — m columns of the 
new Gale dual W will of course be linear combinations of the columns of the old one, but in 
general we will not be able to give an explicit expression for them. We therefore only treat 
some special cases that arise in the context of Ehrhart and chromatic polynomials, and defer 
further discussion to Section 16.21 

6. Applications 

6.1. Ehrhart polynomials. Prom [1], we know that the following inequalities hold for the 
coefficients of iq in the binomial basis: 

fld + ad_i H V Od-s < «o H Vos + Os+i for all < s < [(d - l)/2j . 

For s = 0, the inequality reads < oo + Q^i) and Wd+2 is 

[po - 90, Pi, 0, . . . , 0, -rrf„2)- 

Equation p2|) and Lemma 15.21 thus specialize as follows: 

. ,7r, \{-lY{po-qo)Di,d-{-lfpiroD2.d-To---rd-2 if ^ = 1, 

det Wsn k d\ = { , . J 

^"^ \{-lYpQ---Pk-iDk4-Do,krk-i---rd-2 if 2 < A; < d - 1, 

det W{i^k,i} = (-l)'^9oPo • • • Pk-iDk,iri-i ■ ■ ■ rd-2 - PoDi^kn-i ■ ■ ■ rd-2 {l = dj 

(Here we have used Iverson's notation: |/ = dJ evaluates to 1 if I = d, and to otherwise.) 
Explicit calculation using Lemma 14.11 yields that the coefficient of the leading term r^"^"^ in 
(-l)'^det W^|o,fc,4 is 2(-l)^+i for A; = 1 and (-l)'^+'=+i(d - 2k) for 2 < A; < d - 1. Thus, the 
sign of this coefficient is 

• /r 2d 21/ n TT. ^ f(-l)'^+'=+i foT < k < d/2, 

13 sign (r^'^-^K-l)'' det VFiofedi) = < , 

y > & \y w ) |o,Mi; |(_i)'^+fc ford/2<A;<d. 

We examine the effect that this has on a{J,z). If z G C does not lie in clDo,d;i, the 
first three entries of (jlip already yield two different signs, no matter what sign the last 
determinant takes. Now let z lie inside clDo,d;i, but outside the union of all clDij-i with 
ihj) 7^ i^id)- If {j,k,l} does not contain {0, d}, the ffist three signs of a{J,z) in (fTTj) will 
again contain two different ones. The interesting situation is thus J = {0,k,d}, in which 
case (j{J,z) = |(— 1)'^'+'^+^, (— l)'^signdet VF|o.fc^rf}(2;)}. Combining this with ([13]), we see that 
these signs are different, i.e., z's "last opportunity" J also does not make it an Ehrhart zero, 
if z lies inside the outermost component of the zero locus of det W^|o,fc,d} for < A; < (i/2. 




Figure 6. From top left to bottom right, the vanishing loci of Dq ^ and VFo,fc,d5 
for d = 10, j = 0, I = d, and 1 < k < d — 1. In the first four pictures, the 
outermost oval is only partially shown, but in fact curves around to the right 
and intersects the real axis at a point with large positive coordinate. Thus, 
the points with positive real part just inside P^d;! lie inside the outermost 
component of the zero locus of det M^|q ^ for < k < d/2, but outside all 
components of detM^|o,fc,d} foi" d/2 < k < d. 

6.2. Chromatic polynomials. Let G be a graph on d vertices. The value of the chromatic 
polynomial P{G, t) oi G at z = Iq counts the number of colorings of G with zq colors. The 
chromatic number xiG) is the first positive integer that is not a zero of P{G,t). 

Proposition 6.1. Let G be an undirected graph on d vertices with m edges, k connected com- 
ponents, chromatic number x = x{G), and lo acyclic orientations. Let P{G,z) = 'Ylt=Q^i^i 
be the chromatic polynomial of G expressed in the basis B = {6o, . . . , hd} of P^- 

(a) Let hi = (— so that B is the alternating power basis. Then > for i = 0, . . . ,d, 

= for i = 0, . . . , K — 1, > 0, 0(^-1 = m and a^i = 1. [HI Theorem 2.7]. 

(b) Let hi = Z-, SO that B is the falling factorial basis. Then ai > for i = 0, . . . ,d [HJ 
Theorem 2.1]. 

(c) Let hi = (— SO that B is the alternating rising factorial basis. Then a-i > for 
i = 0, . . . , d m Proposition 2.1]. 

(d) Let hi = (^^^~*) , so that B is the binomial coefficient basis. Then a, > for i = 0, . . . , d, 

Yli=o ~ o,i = ^ precisely for < i < x — 1? o.'nd ad = oj ^ Proposition 4.5], [6]. 

The roots of chromatic polynomials simultaneously satisfy all restrictions implied by these 
nonnegativity conditions. Here we only treat two of these in any detail. 
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6.2.1. The alternating power basis. To evaluate these conditions for the alternating power 
basis, only slight modifications from the power basis case are needed. First, qk = —2x instead 
of Qk = 2x in Proposition 12. H and so 



Z — Z 

with n = k — j — 1. Next, the relations ao = • • • = Ok-i = say that effectively, 



Rk Rk+I ■ ■ ■ Rd 
Ik Ik+1 ■ ■ ■ Id 



( 



\adf 



0, 



so that W starts out with the column (p^i —Qk-.'^'k, 0, . . . , 0)"^. But in the present case of the 
alternating power basis, none of p^, g^, actually depends on k. The matrix W thus stays the 
same, only the effective dimension has dropped io d! = d — k. The discussion in Section 13.21 
still applies, except that the excluded region for roots of -P(G, z) is now the opposite half-open 
sector, i.e., the cone r bounded by the lines of angles ±(1 — 33^)7r. 

We may incorporate the linear equation mad — o-d-i = by appending the row vector 
(0, . . . , 0, — 1, m) of length d — k + 1 to W, and replacing the last two columns of W by their 
linear combination (0, . . . , 0, g, h, m, 1)"^ with g = (m — 2x)(x^ + y^) and h = {m — 2x)2x + 
x'^ + y'^ . The rows of the resulting matrix W represent d' + 1 vectors in IR,"''"^, so any linear 
hyperplane spanned by members of this set is defined by a linear form (pij^k,i- The signs of 
the values of this linear form on the four row vectors W^,w'j,w'f,,W'i are 

<^i,j,k,i = |(-l)'signdetTF'(j_fc^;), (-1)-'+^ signdet l^Jj^fc^,), 

(-l)'^ signdetW;,,^. (-1)^+1 signdet W;,,^. 

(^i,j,k,d-l 



where W(j^ ^.), for instance, is obtained from w' by deleting rows i,j, k. The sets of signs 



|(-l)*signl)j-fc, (-1)^+1 sign (-l)'=signA,i,o| 



tell us that any root allowed by the conditions aj > is also allowed under the additional 
restriction maa = ad~i, so that the set of possible roots does not change under this restriction. 

6.2.2. The binomial coefficient basis. The relations ao = • • • = cl^~i = say that effectively. 



Rx -^X+l • ■ • 
Ix ••• 



0, 



so that W starts out with the column (p^, —q^., r^, 0, . . . , 0)*. The transformation x ^ x -\- x 
maps — g^+j, r^+j) to {pi, —qi,ri), so after this translation the effective dimension has 

dropped to d' = d — x- 

The two afhne linear relations + ■ ■ ■ + a^ = dl and a^ = lo course do not individually 
influence the location of roots, but may be combined to the linear relation Yli=x ~ e'^d = 
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with £ = ^- A Gale dual compatible with this linear relation is the matrix W of size 
(d' + 1) X (d' — 2) with columns 

W = {Vi-Vo, Vd'-S-Vo, XVd'-2- ^Vo) , 

where Vi is the i-th column of W, and the coefficients are A = ed{d — 1) and /i = A — rd'_2- To 
calculate the sets (yi,'j,k,i of signs, we must evaluate the determinant \W\k of the submatrix 
of W obtained by deleting the three rows indexed hy K = k}, say. By multilinearity of 
the determinant, we obtain 

[W]k = Adet - uo, ■ ■ ■ ,Vd'-3 - vq, ^^'-2) - (-1) ~Aidet(fo, ''^i, ...,Ud'-3) 

d'-3 
c=0 

d'-2 

This formula can be evaluated using Lemma 15.21 In Figure [7] we show the zero loci of [W]k 
in the case d = 4 and uj = ^■ 



7. Distribution of random roots 

In closing, we explain a phenomenon encountered several times in the literature [1], [1]: 
The roots of "randomly" generated polynomials with nonnegative coefficients tend to clus- 
ter together in several clumps, and usually lie well inside the region permitted by theory; 
cf. Figure [51 

Our explanation is this: in these simulations, the coefficient vector (oq, . . . ,ad) is usually 
picked uniformly at random from some cube [0, N]'^^^ (except that sometimes the cases ao = 
and ad = are excluded; we will gloss over this minor point). By linearity of expectation, the 
expected value E[f{zo)) of f{zo) = Yli=o aihizo) at a point G C is Y^i=Q E{ai)bi{zo) = 
T Sf=o ^«(-^o)- Thus, as a first approximation, the closer the barycenter /3{zo) = Yli=o ^i(-^o) 
is to zero, i.e., the smaller its absolute value |/3(zo)|, the more likely it is for zq to be a root 
of /! For example, in the case of the binomial coefficient basis, 

by an elementary identity for binomial coefficients. Figure [9] shows the regions where |/3(2;o)| is 
small, together with the roots of several random polynomials. Note that /3(zo) is the Ehrhart 
polynomial of the simplex convjei, . . . , e^, — ei — • • • — e^} by [21 Proposition 1.3]; see also [^. 

Figure [TOl shows the corresponding regions for the rising and falling factorial bases; in the 
case of the power basis (6j = z*), of course (^{zq) = iff 7^ 1 is a d-th root of unity. 

Clearly, the predictive power of this simple model can be easily improved by considering 
additional parameters of the data; however, we will not do this here. 
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Figure 7. From left to right and top to bottom, the zero loci of [VFjii- for 
d = 4 and £ = ^- Here K runs through (^^^'-''^^^ in lexicographic order, except 
that the zero loci coresponding to K = {0, 1, 2} and K = {1, 2, 3} are empty 
and not shown. The last figure combines all the zero loci with the roots of 
500 random polynomials whose coefficients satisfy Yli=o ~ = 0. 
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Figure 9. Left: For d = 6, the contours \P{zo)\ = c for bi = C""^^"') and vary- 
ing c; the innermost contours correspond to the smallest c. Right: additionally, 
the roots of 500 polynomials of degree d whose coefficients with respect to the 
bi are chosen uniformly at random from [0, dl], except that oq, ad ^ 0. 



Appendix 

Proof of Lemma \4^ It suffices to do the calculation for Dn from equation ([5]). So let's expand 
the difference 

n+l n+1 

(14) Yliz - a,){z + a,) - JJ (z - ai)(z + a^), 

1=1 i=l 
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-1 1 2 3 4 5 -6 -5 -4 -3 -2 -1 1 



Figure 10. Contours \P{zo)\ = c for d = 6 and bi = z- (left), hi = (right), 
together with the roots of 1000, respectively 100 random polynomials with 
nonnegative coefficients with respect to these bases. For the rising factorial 
basis, the minima of |/3(2o)| turn out to be real, so they only govern the 
distribution of the real zeros. 



and pick out a term in the expansion with / 'z's and m — I 'z's. The coefficient of this term 
is a sum of terms of the form 



(-1)" 'oii •••ai„_,aji •••aj, - a^, • • • a^^,, (-l)'ai^_,aj, • • • a 



and each of these terms vanishes for ra even. In particular, the term z^^^z^'^^ does not occur, 
which is also easy to see directly. The first nonzero term in (jl4p is then 



2z"z"+i(-ai fln+i) + 2z"+iz"(ai + • • • + a„+i) = 2z"z"(z - z)(ai + • • • + a„+i). 

It is easy to work out ^^^1=1 ^« ~ ^('^ + 1)/^ fo'^ Oj = i + A/2, and this finishes the proof. □ 



Proof of Proposition \4.4[ The curve V = T>n is also described by the equation g = {hi — 
h2)/{z - z), where hi = ]Xi^l{z - ai){z + a,), /12 = nr=i^('^ + "^*)(^ - Oj), and Oj = i + A/2 
with A = d — 2 — n. Thus, D has a singular point if and only if the Jacobi matrix of g 
vanishes at some point of the locus g = 0. Using the chain rule and the relations dz/dx = 1, 
dz/dy = i, dz/dx = 1, dz/dy = —i, we calculate the partial derivatives of g{z) with respect 
to X and y: 

dg{z) - /i2,z + ^1,2 - h2,z 



dx z — z 

dg{z) _ 2^ hi - h2 _^ . - h2,z - hi^z + h2,z 
dy {z — z)'^ z — z 

Here hj^z, hj^z denote the partial derivatives of hj with respect to z, z; by explicit differenti- 
ation, hi^z = Yh=i ^1/ - (^i) and h2,z = Th=i ^2/ (z + a^). 
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To prove that T> has no real singular points, we pick z (z T> and calculate 

% 

dx 



[z] = — ■ 



z — z z — z 

fc, ^/ 1 1 1 1 



z — z ^-^ V z — ~'~ z + a, 2; + Oj z — tti 

1=1 



4^1 E 



n+1 



^ ((x - ai)2 + y2) ((3; + + y2) ' 

For real nonzero z € this expression never vanishes. The same calculation already proves 
the second statement, because a tangent vector to the curve g{x,y) = at a non-singular 
point (xo,yo) is given by ±{-^{xo,yo), |f (xo,yo))- 

We now examine a non-real singular point zq of V. Any such point must satisfy 

= = hi^z{zo) - h2,z{zo)- 

zo - Zq 

The first equation tells us that hi{zo) = h2{zo), so that hi^zizo) — h2^z{zo) = if and only if 
/ii(zo) = (which is incompatible with zq and g{zQ) = 0), or 

n+1 ^ n+1 ^ n+1 

= - - = 2 V 

^ zo - o, ^ 20 + ^ ^0 - 

Writing Zq = xq + zuq and separating the real and imaginary parts in the last expression yields 

n+1 n+1 3 

V - = V ^ = 

^ (xo - a,)2 + y2 _ a,)2 + y2 

But the denominators of these expressions are positive (the Oj and the origin do not lie on V), 
and Ui > for i = 1, . . . , n -|- 1, so we conclude that V has no singular points. □ 
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